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ABSTRACT 

FLAASH  (Fast  Line-of-sight  Atmospheric  Analysis  of  Spectral  Hypercubes)  is  a  first-principles 
atmospheric  correction  algorithm  for  visible  to  shortwave  infrared  (SWIR)  hyperspectral  data.  The 
algorithm  consists  of  two  main  steps.  The  first  is  retrieval  of  atmospheric  parameters,  visibility  (which  is 
related  to  the  aerosol  type  and  distribution)  and  column  water  vapor.  The  second  step  is  solving  the 
radiation  transport  equation  for  the  given  aerosol  and  column  water  and  transformation  to  surface 
reflectmce.  The  focus  of  this  paper  is  on  the  FLAASH  water  vapor  retrieval  algorithm.  Modeled  radiance 
values  in  fire  spectral  region  of  one  water  vapor  absorption  feature  are  calculated  fironi  MODTRAN  4  using 
several  different  water  vapor  amounts  and  are  used  to  generate  a  LoOk-Up  Table  (LUT).  The  water  band 
typically  used  is  1130  nm  but  either  the  940  or  820  nm  band  may  also  be  used.  Measured  radiance  values 
are  compared  to  the  LUT  to  determine  the  column  water  vapor  amount  for  each  pixel  in  the  scene.  We 
coirpare  the  results  of  water  retrievals  for  each  of  these  bands  and  also  the  results  of  flieir  corresponding 
reflectance  retrievals. 

Keywords:  hyperspectral,  atmospheric  correction,  water  vapor  retrieval,  FLAASH 

The  Fast  Line-of-sight  Atmospheric  Analysis  of  Spectral  Hypercubes  (FLAASH)  atmospheric 
correction  algorithm/code  is  a  software  package  developed  by  the  Air  Force  Research  Laboratory,  Space 
Vehicles  Directorate  (AFRLA^S),  Hanscom  AFB  and  Spectral  Sciences,  Inc.  (SSI)  to  support  the  analyses 
of  visible  to  shortwave  infrared  (Vis  -  SWIR)  hyperspectral  and  multispectral  imaging  sensors.  The 
algorithm  derives  its  first-principles  physics-based  calculations  from  the  MODTRAN4  (Berk  et  al,  1996; 
Anderson  et  al,  2000)  radiative  transfer  code.  The  main  objective  of  FLAASH  is  to  eliminate  atmospheric 
effects  caused  by  molecular  and  particulate  scattering  and  absorption  from  the  ‘radiance-at-detector’ 
measurements  in  order  to  retrieve  ‘reflectance-at-surface’  values.  In  order  to  achieve  this  objective, 
FLAASH  enploys  two  main  steps.  The  first  is  retrieval  of  atmospheric  parameters,  visibility  (which  is 
related  to  the  aerosol  type  and  distribution)  and  column  water  vapor.  Since  current  methods  allow  aerosol 
retrieval  over  a  very  limited  set  of  surface  types  (dark  land  pixels),  only  an  average  visibility  is  obtained  for 
a  scene.  On  the  other  hand,  the  spectral  signature  of  water  vapor  is  sufficiently  distinct  that  the  column 
amount  may  be  retrieved  on  a  pixel-by-pixel  basis.  The  second  step  is  solving  the  radiation  transport 
equation  for  the  given  aerosol  and  column  water  and  transformation  to  surface  reflectance.  The  emphasis  of 
this  paper  is  on  the  application  of  the  FLAASH  water  vapor  determination  algorithm  and  subsequent 
retrieval  of  surface  reflectance  spectra  using  AVIRIS  hj^erspectral  radiance  measurements  as  the  input. 

*Gerald.Felde@hanscom.af.mil;  phone  781-377-3136;  fax  781-377-3138 


2.  DESCMPnON  OF  FLAASH 


MODTRAN  calcukted  qjectral  radiance  outpnte  are  central  to  FLAASH  operation.  The  MODTRAN 
calcuktions  use  die  viewing  and  solar  angles  and  the  mean  surfece  elevation  of  the  measured  scene.  In 
addition^  selection  of  reasonable  MODTRAN  standard  atmosphere  and  aerosol  models  should  be  made. 
The  outputs  from  MODTRAN  are  convolved  with  the  sensor  filter  functions  so  tiiat  flie  radiance 
simulations  have  the  same  spectral  characteristics  (i.e.,  clmnnel  center  wavelengths  and  FWHMs)  as  the 
measurements.  A  typical  FLAASH  run  contains  three  sets  of  MODTRAN  calculations.  The  jSrst  set  is 
required  for  die  preliminary  column  water  retrieval  of  each  pixel  in  the  scene.  A  i^er  estimate  of  the  scene 
visibility  is  an  additional  input  needed  for  this  set.  The  mean  value  of  the  retrieved  column  water  is  used  as 
input  to  the  second  set,  which  is  required  for  the  scene  average  visibility  retrieval.  TTiis  vkibility  value  is 
used  as  input  to  the  third  set,  which  is  used  for  the  final  column  waiter  vapor  deterrhination  and  reflectance 
spectra  retrieval  of  each  pixel  in  the  scene.  For  the  water  column  md  visibility  MbDTRAN  calculations, 
the  parameter  of  interest  (for  visibility  it  is  actually  inverse  visibility)  is  varied  in  even  steps  over  the  likely 
range  of  values  in  the  cube.  The  p'ound  albedo  is  also  varied  over  die  values  of  0,  0.5,  and  1,  so  there  are 
two  nested  loops.  The  spectral  range  of  each  set  of  MODTRAN  calculations  must  contain  die  wavelengths 
required  for  the  task  at  hand  -  the  relevant  water  bands  for  die  preliminary  water  column  determination,  die 
bands  inquired  for  die  visibility  determination,  and  the  full  range  of  the  instrument  for  die  final  water 
column  determination  and  reflectance  spectra  retrieval. 

FLAASH  processes  mdiance  images  with  Rectal  coverage  fi-om  wavelengths  in  the  solar  reflective 
regime  where  themml  emission  can  be  neglected.  For  this  situation  the  spectral  radiance  L*  at  a  sensor 
pixel  imy  be  parameterrod  as  (Vermote  et  a/.,  1994,  Staenz  et  al^  1996,  Adler-Golden  et  aL,  1998, 
Matthew  ef  a/.,  2000) 

L*=Ap/(l-PeS)+Bpe/(LpeS)  +  L^  (1) 

where  p  is  the  pixel  surface  reflectance,  pe  is  an  average  surface  reflectance  for  the  surrounding  region,  S  is 
the  spherical  albedo  of  the  atmosphere  (capturing  the  backscattered  surface-reflected  photons),  L*^  is  the 
radiance  bacfecattered  by  die  atmosphere  without  reaching  the  surface,  and  A  and  B  are  surface 
independent  coefficients  diat  vary  with  atmospheric  and  geometric  conditions.  The  spatially  averaged 
reflectance,  Pe,  is  used  to  account  for  **adjacency  effects**,  i.e.,  radiance  contributions  dmt,  because  of 
atmospheric  scattering,  originate  firom  parts  of  the  surface  hot  in  the  direct  line  of  sight  between  sensor  and 
targeted  pixel.  All  of  die  variables  are  implicidy  wavelen^  dependent.  The  first  term  in  Equation  (1) 
corresponds  to  the  radiance  reaching  the  surface  (from  bodi  sky-shine  and  direct  solar  illumination)  that  is 
backscattered  directly  into  the  sensor,  while  the  second  term  corresponds  to  the  radiance  from  die  surface 
that  is  re-scattered  by  the  atmosphere  into  the  sensor. 

After  the  atmospheric  retrievals  of  visibility  and  colunm  water  vapor  (discussed  later)  are 
performed.  Equation  (1)  is  solved  for  the  pixel  surface  reflectances  in  all  of  foe  sensor  channels.  Ihe 
solution  is  based  on  a  method  in  which  foe  spatially  averaged  radiance  image  L*c  is  used  to  estimate  foe 
averaged  reflectance  pg.  This  is  done  using  an  approxumte  equation  derived  from  Eqimtion  (1): 

L%  =  (A+B)pg/(l-peS)  +  L%  (2) 

As  in  Equation  (1),  foe  variables  are  all  wavelength  (or  equivalently,  channel)  dependent,  L*e  is  obtained 
from  foe  measured  radiance  values,  L*,  by  spatial  averaging  of  foe  data  using  a  point-spread  function  that 
describes  foe  atmospheric  scattering  of  ground  reflected  photons  into  foe  sensor  pixel.  The  ^atial 
averaging  excludes  foe  measured  values  for  cloud-filled  pixels  because  clouds  do  not  contribute  to  foe  near¬ 
surface  pixel  cross-talk.  These  pixels  are  rqilaced  by  foe  scene  average  radiance  of  all  cloud-free  pixels. 
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The  cloudy  pixels  are  found  using  a  combination  of  brightness,  color  balance,  and  column  water  vapor  tests 
(Ackerman  et  ah^  1998  and  Borel  et  ah,  1999).  Starting  from  the  known  L*e  for  each  pixel  and  the 
coefficients  A,  B,  S  and  L*a  (which  are  extracted  from  the  MODTRAN  radiance  simulations).  Equation  (2) 
is  solved  for  Pe.  The  result  is  inserted  into  Equation  (1),  which  is  then  solved  for  the  reflectance  p. 

There  is  a  final  optional  step  available  in  FLAASH  that  can  be  applied  to  the  retrieved  reflectance 
cube.  It  is  called  spectral  polishing,  a  term  coined  by  Boardman  (1998)  to  describe  a  mathematical 
renormalization  method  for  removing  artifacts  from  reflectance  spectra  using  only  the  data  itself.  When 
properly  in^lemented,  poUshing  dramatically  reduces  spurious,  systematic  spectral  stracture  due  to 
wavelength  registration  errors  and  molecular  absorption  residuals  while  leaving  true  spectral  features  of  the 
suiface  intact.  It  is  important  to  note  that  FLAASH  has  an  automated  wavelen^  recalibration  algorithm, 
which  can  also  be  Used  to  improve  the  r  reflectaUce  spectra.  This  algorithm  uses  atmospheric 

molecular  absorption  features  (e  g.,  oxygen,  carbon  dioxide,  water  vapor)  to  quantify  wavelength  error  by 
fitting  an  absoiption  band  shape  in  the  measured  radiance  spectra  to  MODTRAN  model  calculations  (Felde 
et  a/.  ,  2003).  If  the  wayelength  recalibration  option  is  selected,  it  is  implemented  in  ■  the  FLAASH 
processing  flow  immediately  after  the  preliminaEry  water  vapor  retrieval.  ;  ^  ^  ^ 

The  visibility  or  aerosol  overburden  (optical  column,  not  vertical  structure)  is  currently  addressed 
through  a  general  reflectance  ratio-based  algorithm,  retrieving  only  the  scene-average  aerosol  amount.  The 
algorithm  is  based  on  the  enpirical  observation  by  Kaufinan  et  al  (1997)  that  for  natural  dark  terrain 
(primarily  vegetation)  the  reflectances  at  certain  wavelengths,  such  as  660  nm  and  2100  nm,  are  in  nearly  a 
fixed  ratio  (for  these  wavelengths,  the  ratio  is  --0.45).  The  Equation  (1)  and  (2)  calculations  are  iterated 
over  a  fange  of  yi$ibility  values  using  the  appropriate  Set  of  MOD^^  radiance  simulations  (discussed 
above).  For  each  Aasibilify/^  scene-average  660  mn  aud  2100  niti  reflectances  for  the  dark  pixels  are 
retrieved,  and  the  “best”  esiimatei  of  visibility  is  interpolated  by  matching  the  ratio  to  6.45. 

3.  GOLtJMN^^^ 

Water  vapor  provides  the  strongest  molecular  absorption  featines  in  the  visible  to  SWIR  range.  In  fact, 
water  vapor  absorption  is  so  strong  in  the  vicinity  of  1400  and  1900  nm  that  the  atmosphere  is  opaque  and 
thus  it  is  physically  iirpossible  to  retrieve  water  vapor  and  surface  reflectance  in  these  wavelengfli  regions. 
Water  vapor  absoiption  is  also  centered  near  820,  940,  and  l  lSO  nm.  However,  at  these  wavelengths,  the 
atrriosphere  is  t>^ically  not  opaque  and  so  water  and  surface  reflect&ce  retrieval  is  possible^ 

The  physical  basis  of  the  FLAASH  water  vapor  retrieval  algorithm  is  illustrated  in  Fig.  1.  This  figure 
shows  the  effect  Of  varyhig  iauhoimts^^^^p^^  water  vapor  (0,  2006,  4000,  ...  TOOOO  on  the 

radiance  spectra  in  the  vicinity  of  the  820,  940,  and  1130  nm  water  absorption  regions.  The  absorption 
“valleys”  get  deeper  with  increasing  amounts  of  water  yapor.  These  radiance  spectra  were  generated  from  5 
cm“*  resolution  MODTRAN  calculations,  which  were  convolved  to  correspond  to  the  AVIRIS  channel 
center  wavelengths  and  FWHMs.  The  AVIRIS  channel  center  wavelengths  are  -10  nm  apart  and  the 
FWHMs  are  also  -lOnm.  The  MODTRAN  calculations  used  the  viewing  and  solar  conditions  (56  degree 
solar  zenith  angle)  corresponding  to  an  AVIRIS  radiance  data  cube  for  Fort  A.P.  Hill,  VA,  which  was 
analyzed  with  FLAASH  (this  analysis  is  discussed  later).  Also,  the  calculations  used  a  surface  albedo  of 
0.5,  the  FLAASH  retrieved  scene-average  visibility  of  33  km,  the  rural  aerosol  model,  the  Isaacs  2-stream 
multiscatter  model,  and  the  tropical  atmospheric  profiles.  Eleven  sets  of  MODTRAN  calculations  were 
done,  keeping  all  parameters  constant,  except  for  scaling  the  tropical  moisture  profile  to  correspond  to 
column  water  values  of  0,  1000,  2000,  ...  10000  atm-cm.  (It  should  be  noted  that  the  0  atm-cm  simulation 
actually  was  based  on  a  scale  factor  of  0.01  and  thus  in  reality  used  a  moisture  profile  containing  51  atm- 
cm  column  water  vapor.  This  is  why  the  curve  labeled  0  water  in  Fig.  1  has  small  valleys  near  820,  940, 
and  1130  nm.  The  curve  for  an  exact  value  of  0  water  would  not  have  these  small  valleys.)  For  clarity,  only 
six  of  the  fesulting  eleven  isimulated  radiance  spectra  were  plotted  in  ^  F^^  I.  There  are  six  model 
atmospheres  available  in  MODTRAN  (Table  1)  and  their  column  water  vapor  values  range  from  -  500  to 
5000  atm-cm.  In  the  Fig.  1  MODTRAN  radiance  simulations,  the  maximum  water  value  used  was  twice 


that  of  tbe  typical  tropical  moisture  profile  available  in  MODTRAN.  TTiis  was  done  to  account  for 
extremely  mokt  atmospheres,  which  can  occur  in  nature. 


775  825  875  925  975  1025  1075  1125  1175  1225 

wavelength  (nm) 

Fig.  1 .  MODTOAN  simulated  radiance  s|^ctta  for  varying  amounts  of  column  water  vapor. 


Table  1.  Column  Water  Vapor  Values  of  MODTRAN’s  Atmospheric  Models. 


Atmospheric  Model 

Column  Water  Vapor 
(atm-cm) 

sub-arctic  winter 

518 

mid-latitude  winter 

1060 

U.S.  standard 

1762 

sub-arctic  summer 

2589 

mid-latitude  summer 

3636 

tropical 

5119 

In  order  to  quantify  and  conq>are  the  sensitivity  to  water  vapor  in  the  820,  940,  and  1130  nm  regions 
using  the  simulated  radiance  spectra,  the  depth  of  each  absorption  valley  for  each  column  water  vapor 
amount  was  confuted.  This  was  done  by  dividing  each  absorption  band  minimum  radiance  value  by  the 
zero  water  radiance  value  at  the  same  wavelength.  The  resulting  values  were  multiplied  by  100  in  order  to 
tiansform  them  into  percentages.  The  final  resulte  are  plotted  in  Fig.  2.  There  are  three  curves  of  dqsth  of 
absorption  valley  vs.  column  water  vapor  -  one  for  each  of  the  three  absorption  regions.  It  should  be  noted 
that  the  greater  Ae  depth,  the  smaller  the  value  on  the  vertical  axis.  For  the  column  water  vapor  range  of  0 
to  ~  4000  atm-cm,  the  940  and  1 130  imn  curves  are  similar.  Both  show  a  rapid  increase  in  valley  depfii  with 
increasing  water  vapor,  which  indicates  a  strong  sensitivity  to  water  vapor.  On  the  odier  hand,  for  this  same 
water  range,  the  820  nm  curve  shows  a  much  more  gradual  increase  in  valley  depth  with  increasing  water, 
which  indicates  a  much  weaker  sensitivity  to  water.  For  water  vapor  values  greater  than  4000  atm-cm,  die 
magnitude  of  the  slopes  of  all  three  curves  is  small  conq^ared  to  their  respective  slopes  for  water  values  less 
than  4000  atm-cm.  This  indicates  less  sensitivity  to  water  vapor  in  all  diree  absorption  regions  for  water 


vapor  values  greater  than  4000  atm-cm.  The  difference  between  the  4000  and  7000  atm-cm  percent 
radiance  values  are  9.10,  8.40,  and  3.90  for  the  820,  940,  and  1130  nm  ciirves,  respectively.  This  indicates 
similar  sensitivity  to  water  vapor  at  820  and  940  nm  for  this  range  of  vapor  values  but  weaker  sensitivity  at 
1130  nm.  The  difference  between  the  7000  and  10000  atm-cm  percent  radiance  values  are  6.48,  3.41,  and 
0.87  for  the  820,  940,  and  1130  nm  curves,  respectively.  Thus  the  greatest  sensitivity  to  water  vapor  for  this 
water  range  is  at  820  nm.  At  1130  nm,  the  radiance  values  are  nearly  zero  for  this  water  range,  which 
indicates  that  the  atmosphere  is  essentially  opaque  and  thus  has  no  sensitivity  to  water  vapor.  The  above 
discussion  points  out  that  either  the  1 130  or  940  nm  water  vapor  band  is  the  best  choice  for  water  retrieval 
for  non-tropical  moisture  conditions  (recall  that  the  tiopical  moisture  profile  in  MODTRAN  has  a  column 
water  vapor  value  of  51 19  atm-cm,  see  Table  1)  and  that  either  the  940  or  820  nm  band  is  the  best  choice 
for  retrieval  under  very  moist  conditions. 


colunfin  water  vapor  (atm-cm) 

Fig.  2.  Absorption  yalley  depths  for  three  water  vapor  bands  as  a  function  of  column  water  vapor  amount. 

ITie  principal  component  t>f  the  FLAASH  column  water  Vapor  retiieval  algorithm  is  the  strong 
dependence/correlation  of  water  vapor  column  amount  to  the  ratio  of  ‘reference’  (the  shoulders  of  the  water 
absorption  band)  and  ‘absorption’  (center  of  the  same  absorption  band)  radiances.  The  retrieved  water 
vapor  also  depends  somewhat  on  the  absolute  values  of  these  radiances,  which  vary  most  directly  with  the 
surface  reflectance.  This  dependence  arises  because  the  amount  of  water  absorption  in  the  atmospherically- 
scattered  and  surface-reflected  radiance  components  is  slightly  different;  the  absorption  is  generally  smaller 
for  the  atmospherically  scattered  photons,  which  avoid  the  high  concentration  of  water  vapor  close  to  the 
ground.  A  Look-Up  Table  (LUT)  is  then  built  around  these  dependencies. 

The  FLAASH  column  water  vapor  retrieval  algorithm  uses  Equation  (1)  but  for  simplicity  ignores  the 
adjacency  effect,  replacing  Pe  with  p.  This  is  a  good  approximation  because  the  wavelengths  required  for 
water  retrieval  are  fairly  long  and  thus  atmospheric  scattering  is  minimal.  The  MODTRAN  generated  A, 
B,  S,  and  L*a  coefficients  extracted  for  each  channel  and  column  water  amount  are  used  to  simulate  a  set  of 
absorption  and  reference  radiances,  L*.  This  is  done  for  an  even  grid  of  surface  reflectance  values  from  0  to 
1,  inclusive,  spaced  0.01  apart.  The  ratio  of  the  reference  to  absorption  bands  is  computed  for  each  grid 
point.  The  result  is  a  set  of  triples,  consisting  of  {column  water  vapor,  reference  rachaiice,  ratio}  that  span 
both  the  entire  expected  range  of  water  column  and  surface  reflectances.  These  values  are  initially  evenly 


spaced  in  the  dependent  variable  to  be  retrieved  (water  vapor)  but  not  in  the  independent  variables 
(measured  radiance  -  reference  and  ratio). 

The  triples  are  transformed  into  an  even  grid  in  reference  radiance  md  ratio.  The  msulting  2« 
dimensional  LUT  has  reference  radiance  and  ratio  as  independent  variables  and  the  water  vapor  as  the 
dependent  variable.  This  LUT  allow  a  direct  table  look-iq>  of  water  column^  given  measured  reference 
radiance  and  ratio,  and  thus  it  is  searched  to  retrieve  the  column  water  vapor  value  for  each  pixel.  Jn 
occasional  instances,  die  reference  radiance  or  ratio  value  may  lie  outside  the  LUT.  In  such  cases  the  pixel 
is  flagged  to  denote  LUT  tailum,  and  the  water  vapor  for  that  pixel  is  then  assigned  the  average  of  the  valid 
retrievals  for  the  scene.  The  water  determination  is  a  necessary  precursor  to  final  reflectance  retrieval  of 
each  individual  pixel  because  it  permits  calculation  of  the  correct  values  for  the  coefficients  in  Equations 
(2)and(l). 

4.  FI^SH  WATER  Al^  REFLECTANCE  RET^ 

Three  cloud-free  AVIRIS  radiance  ckta  cubes  collected  under  substanrially  different  atmospheric 
moisture  conditions  ware  selected  for  analysis  with  FLAASH:  (1)  Fort  Huachuca,  AZ  on  12  November 
1999,  (2)  Fort  A.P,  Hill,  VA  on  8  Noveniber  1999,  and  (3)  Can^  Ripley,  MN  on  10  August  2000.  For  each 
radiance  data  cube,  three  different  column  water  vapor  retrievals  were  performed,  using  the  1 130, 940,  and 
820  nm  absorption  features,  respectively.  Also,  three  different  reflectance  cubes  were  retrieved  from  each 
radiance  cube  using  the  1 130, 940,  and  820  nm  water  resulte,  respectively.  The  following  FLAASH  options 
were  enqsloyed:  autonmted  wavelength  ^calibration,  adjacency  correction,  and  spectral  polishing. 

The  water  vapor  column  retrieval  results  for  these  three  cubes  are  given  in  Table  2.  Invalid  pixels  are 
those  whose  water  reference  radiance  or  ratio  (of  reference  radiance  to  absorption  radiance)  value  lie 
outside  the  LUT  and  thus  retrieval  of  water  vapor  was  not  possible.  The  majority  of  invalid  pixels  are 
typically  those  containing  surface  water  (elg.,  lakes).  Often  it  is  not  possible  to  determine  water  vapor  over 
surface  water  because  of  its  near  zero  radiance  values  at  the  vravelengths  corresponding  to  the  vrater  vapor 
absorption  features,  which  is  a  lesult  of  the  very  small  reflectance  of  liquid  water  in  this  spectral  legion. 
The  Fort  Huachuca  and  Camp  Ripley  scenes  both  contain  no  invalid  pixels,  which  is  consistent  with  the 
fact  that  there  is  no  surface  water  in  these  scenes.  On  the  other  hand,  the  ~  1  %  invalid  pixels  for  the  Fort 
A.P.  Hill  cube  corresponds  to  it  having  -- 1  %  surface  water  pixels  in  the  scene.  The  average  retrieved 
column  water  vapor  value  obtained  from  using  a  particular  band  for  Canq)  Ripley  is  2  times  that  for  Fort 
A.P.  Hill,  and  the  water  value  for  each  band  for  Fort  A.P.  HiU  is  ~  2.5  times  that  for  Fort  Huachuca.  For 
each  scene,  the  average  retrieved  water  vapor  values  are  similar  for  the  1 130  and  940  nm  absorption  bands, 
but  are  considerably  larger  for  the  820  nm  band.  Ihis  is  consistent  with  the  observations  obtained  from 
Figs,  1  and  2  that  for  atmospheres  containing  less  than  -  4000  atm-cm  column  water  vapor,  the  water  vapor 
column  sensitivity  is  much  stronger  for  the  1130  and  940  nm  absorption  bands  con^ared  to  the  semitivity 
for  file  820  nm  band.  Also,  this  suggesta  tlmt  die  retrieved  water  values  obtained  from  eifrter  the  1 130  or 
940  nm  region  should  be  much  closer  to  die  trae  value  than  the  retrieved  value  obtained  from  the  820  nm 
region.  This  hypothesis  was  verified  by  comparison  of  individual  pixel  reflectance  spectra  for  each  scene 
from  the  three  Afferent  reflectance  cubes,  which  were  retrieved  from  each  radiance  cube  using  the  1130, 
940,  and  820  nm  water  resulte,  respectively.  The  comparisons  showed  that  some  residual  water  vapor 
signal  ws  evident  in  the  reflectance  spectra  obtained  with  the  820  nm  water  results  but  not  in  the 
reflectance  spectra  obtained  with  die  other  two  ban^. 

Examples  of  some  reflectance  spectra  conqiarisons  for  two  pixels  from  the  Fort  A.P.  Hill  scene  are 
given  in  Figs.  3  and  4.  Fig.  3  is  a  building  pixel  and  Fig.  4  is  a  grass  pixel.  The  retrieved  column  water 
vapor  amounte  for  each  of  these  individual  pixels  were  almost  exactly  the  same  as  the  scene  average  values 
listed  in  Table  2.  The  reflectance  spectra  obtained  using  the  1130  and  940  nm  regiom  for  water 
determination  for  each  of  these  two  pixels  was  essentially  identical  because'  the  water  values  were  nearly 
the  same.  Thus,  a  single  curve  is  used  in  each  figure  for  these  two  spectra.  The  FLAASH  reflectance 
spectra  using  die  1 130  or  940  nm  water  vapor  retrievals  are  good,  i.e.,  no  spikes  or  valleys  at  820,  940,  or 


1130  nm.  On  the  other  hand,  the  reflectance  spectra  using  the  820  lun  water  vapor  retrievals  contain 
noticeable  over-correction  (i.e.,  spikes)  at  940  and  1130  nm,  and  also  on  the  wings  of  the  two  opaque 
absorption  regions,  which  are  centered  near  1400  and  1900  nm.  Thir  indicates  that  the  amount  of  water 
retrieved  using  the  820  nm  absorption  feature  was  too  large.  It  should  be  noted  that  there  is  no  spike  on  this 
curve  at  820  nm.  The  820  nm  absorption  feature  is  considerably  weaker  than  either  the  940  or  1 130  feature 
(see  Fig.  1).  Apparently,  the  17  %  larger  water  vapor  value  obtained  from  the  820  nm  absorption  isn’t 
enough  of  a  difference  to  adversely  affect  the  the  reflectance  retrieval  in  the  820  nm  region. 


Table  2.  FLAASH  Water  Vapor  Column  Retrieval  Results. 


Location 

Water  Vapor 
Absorption 
FeatureUsed 

Standard 

Deviation 

Percent  Invalid 
Pixels 

Ft.  Huachuca 

1130  nin 

588 

0 

940imi 

573 

44 

0 

820  nm 

711 

35 

0 

Ft.  A.P.  Hill 

1130nm 

1516 

136 

1.06 

940  nm 

1505 

179 

1.18 

820  nm 

1763 

147 

1.07 

Camp  Ripley 

1130  nm 

2904 

177 

0 

940  nm 

3206 

185 

0 

820  nm 

3849 

111 

0 

Fig.  3.  Building  reflectance  spectra  retrieved  by  FLAASH  using  the  1130  (or  940)  and  820  nm  regions  for  water  vapor 
column  determination. 


Fig*  4.  Grass  reflectance  spectra  retrieved  by  FLAASH  using  the  1 130  (or  940)  and  820  nm  regions  for  water  vapor 
column  determination* 

In  order  to  help  gain  a  better  understanding  of  the  effect  of  water  vapor  column  anK)unt  on  the  FLAASH 
retrieved  reflectance  spectra,  retrievals  were  done  for  die  Fort  A*P.  Hill  cube  using  three  vmter  amounts: 
O.lx,  0.4x,  and  0.7x  die  MODTRAN  mid-latitude  summer  profile,  which  contains  3636  atm-cm  of  water 
vapor.  Reflectance  spectra  resulte  me  shown  in  Figs.  5  and  6  for  the  same  two  pixels  used  for  Figs.  3  and  4, 
respectively  The  0.4  scale  fector  yields  a  column  water  value  approximately  the  same  as  was  correctly 
retrieved  using  either  the  1 130  or  940  nm  absoiption  band.  Thus,  the  reference  (i.e.,  coirect)  spectra  in  Figs. 
5  and  6  k  the  one  obtained  with  the  0.4  water  vapor  scale  factor.  The  0.1  and  0.7  scale  factors  give  column 
values  that  are  75  %  smaller  and  larger  compared  to  die  correct  valuer  respectively.  The  reflectaace  spectra 
in  Figs.  Sand  6  obtained  using  the  0.1  water  vapor  scale  factor  contain  large  valleys  (under-correction 
due  to  use  of  a  too  small  water  value)  at  940  and  1130nm;as  well  as,  on  the  wings  of  the  two  opaque 
absoiption  regions,  which  are  centered  near  1400  and  1900  nmu  On  the  other  hand,  spectra  obtained  using 
the  0,7  scale  factor  contain  large  spikes  (over-correction  due  to  use  of  a  too  large  water  value)  at  die  same 
wavelength  locations  of  the  valleys  of  the  other  curve.  Each  of  the  two  curves  containing  peaks  has  a 
larger  one  at  1 130  nm  compared  to  940  nm.  This  is  also  true  of  each  of  die  two  curves  containing  valleys.  It 
should  be  noted  that  the  two  curves  containing  peaks  don’t  have  any  at  820  nm  but  the  curves  are  a  bit 
higher  than  the  reference  spectra  in  this  wavelength  region.  Also,  the  two  curves  containing  valleys  don’t 
have  any  at  820  nm  but  the  curves  am  a  little  lower  than  die  reference  spectra  in  this  spectral  region.  These 
observations  show  that  if  the  column  water  vapor  value  used  for  the  FLAASH  reflectance  retrieval  is 
incorrect,  then  die  effect  is  larger  at  1 130  nm  than  at  940  nm,  while  the  effect  is  minimal  at  820  nm.  This  is 
consistent  with  the  observation  that  the  1130  nm  absorption  feature  is  more  sensitive  to  water  vapor  than 
the  940  nm  feature  and  that  the  820  nm  band  only  has  small  sensitivity  to  water  when  the  true  column 
water  vapor  value  is  -  1500  nm  (see  Figs.  1  and  2). 
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Building  reflectance  spectra  retrieved  by  FLAASH  using  three  very  different  water 
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5.  CONCLUSIONS 


FLAASH  retrieved  column  wato  vapor  amounte  for  tibree  AVTRIS  radiance  cubes,  which  were 
measured  under  low  to  moderate  moisture  conditiom  (<  -  4000  atm-cm),  are  similar  using  the  1 130  or  940 
nm  absorption  region,  but  are  20  %  larger  using  the  820  nm  region.  The  FLAASH  reflectance  spectra 
results  for  fliese  data  cubes  ming  either  the  1 130  or  940  nm  water  retrievals  are  good,  but  over-correction  in 
the  940  and  1130  nm  water  vapor  absorption  regions  occurs  using  die  820  nm  water  retrievals:  These 
results  are  consistent  with  information  obtained  from  a  column  water  vapor  semitivity  analysis  of 
MODTRAN  calculated  radiance  values,  which  shows  tlmt  the  820  nm  band  im  limited  sei^itivity,  while 
the  940  and  1 130  nm  have  hi^  sensitivity  for  low  to  moderate  amotmts  of  water  vapor.  This  analysis  also 
shows  for  high  values  of  water  vapor  diat  the  sensitivity  is  best  at  820  nm;  however,  it  is  quite  low  in  all 
three  absorption  regions.  There  is  a  need  to  examine  an  expanded  set  of  measured  radiance  cubes,  which 
cover  the  con^Iete  range  of  possible  moisture  conditions,  to  better  quantify  the  differences  in  FLAASH 
water  vapor  column  retrievals  using  the  820,  940,  and  1130  nm  regions  and  also  the  differences  in  the 
corresponding  FLAASH  retrieved  reflectance  ^ectra. 
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